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^ ' The spectral density of the D-meson in the nuclear environment is studied 

within a self-consistent coupled-channel approach assuming a separable po- 
tential for the bare meson-baryon interaction. The DN interaction, described 
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, through a G-matrix, generates dynamically the Ac (2593) resonance. This 

resonance is the charm counterpart of the A (1405) resonance generated from 
the s-wave KN interaction in the 1=0 channel. The medium modification 
j. ■ of the D-meson spectral density due to the Pauli blocking of intermediate 

^ ■ states as well as due to the dressing of the D-mesons, nucleons and pions is 



o 
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investigated. We observe that the inclusion of coupled-channel effects and the 



! self-consistent dressing of the D-meson results in an overall reduction of the 



in-medium D-meson changes compared to previous work which neglect those 
effects. 
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^ ■ I. INTRODUCTION 



The study of the properties of hadrons in a hot and dense medium is an important 
problem in strong interaction physics. It is a topic of active research interest as it has direct 
imphcations in heavy-ion collision experiments, as well as in the study of astrophysical 
compact objects like neutron stars. The medium modifications of hadrons are probed in the 
relativistic heavy-ion collision experiments. In particular, at CERN-SPS, the experimentally 
observed dilepton spectra [1,2] have been attributed to the medium modifications of the 
spectral properties of the vector mesons, especially of the p-meson [3-7] and can not be 
explained by the vacuum hadronic properties. Furthermore, the production from nuclear 
collisions at GSI-SIS energies of 1-2 AGeV have shown that in-medium properties of the 
kaons have been seen in the collective flow pattern of mesons as well as in the abundancy 
and spectra of antikaons [4,8-15]. The medium modifications of D{D) mesons, which show 
analogy to the K{K) mesons resulting from replacing s-quark (s-antiquark) by the c-quark 
(c-antiquark), have also become a subject of recent interest [16-21]. The medium changes 
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for D-mesons can have important consequences for the open charm enhancement in nucleus- 
nucleus collisions [22] as well as for J/^ suppression as observed at the SPS [23]. The NA50 
Collaboration has claimed to see an open charm enhancement by up to a factor of three in 
central Pb + Pb coUisions at 158 A-GeV [24]. The medium modifications of the open charm 
mesons {D{D)) can also modify the high mass (M > 2 GeV) dilepton spectra [25] since 
they can be produced abundantly in high energy heavy-ion collisions. Transverse momentum 
spectra of electrons from Au + Au collisions at ^/s = 130 GeV have been measured at 
midrapidity by the PHENIX experiment at RHIC [26]. The spectra show an excess above 
the background from decays of hght hadrons and photon conversion. The observed signal is 
consistent with that expected from semi-leptonic decays of charmed mesons. 

In high energy heavy- ion collisions at RHIC (^/s ~ 200 GeV), an appreciable contribution 
of J/\E' suppression is expected to be due to the formation of a quark-gluon Plasma (QGP) 
[27]. However, the medium modification of D-mesons should modify the J/"^ absorption in 
hot and dense nuclear medium and can provide a possible explanation for J/^ suppression. 
The effect of the hadron absorption of J/^'s is found to be not neghgible [28-30]. It is thus 
of importance to understand the interactions of the D-mesons in the hadronic medium for 
J/^ production as well. 

At finite densities, the medium modification of the D-meson mass has been studied 
using the QCD sum rule (QSR) approach [16,19]. Due to the presence of a light quark in 
the D-meson, the mass modification of the £)-meson has a large contribution from the light 
quark condensates. The large shift of mass (~ 50 MeV at nuclear matter density, po) in the 
Z^-mesons is originated from the contribution of the mc{qq)N term in the operator product 
expansion. The quark-meson coupling (QMC) model predicts a mass drop of the D-meson 
to be of the order of 60 MeV at p = po [21], which is very similar to the value obtained 
in the QCD sum rule calculations of Refs. [16,19]. At finite temperatures, the studies of 
quarkonium dissociation [20,31] using heavy quark potential from lattice QCD [32] suggest 
a similar drop of the D-meson mass. 

In a recent work, the mass modification of the D-meson in hot and dense matter, arising 
due to its interaction with the light hadron sector, was studied in a chiral effective model 
[33]. The chiral SU(3) model, used to study the hadronic properties in the hot hyperonic 
matter [34] was generalized to SU(4) to include the charmed mesons. Interactions of the 
D-meson with the light hadron sector were derived to study the in-medium mass of the 
D-meson. 

In all these investigations, the spectral features of a D-meson embedded in hot and dense 
matter have not been studied. Coupled-channel effects as well as dressing of intermediate 
propagators have been completely ignored, which turned out to be crucial for describing the 
strange counter part of the D-meson, the K meson in the nuclear medium. In the present in- 
vestigation, the spectral density for the D-meson is calculated using a self-consistent coupled- 
channel G-matrix calculation. The medium effects like the Pauli blocking on the nucleons 
in the intermediate states and the effects from the dressing of the D-mesons, nucleons as 
well as pions are investigated. The coupled-channel formalism generates dynamically the 
Ac(2593) resonance, like the A(1405) resonance for the case of antikaons. 

We have organized the present paper as follows. In Sect. II we first review our formalism. 
Our results are presented and discussed in Sect. HI and the concluding remarks are given 
in Sect. IV. 
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II. FORMALISM 



In this section we present the formahsm to obtain the self-energy or single-particle poten- 
tial of a D-meson embedded in infinite symmetric nuclear matter. This self-energy accounts 
for the interaction of the D-meson with the nucleons and its calculation requires the knowl- 
edge of the in-medium DN interaction, which will be described by a G-matrix. The medium 

effects incorporated in this G-matrix include the Pauli blocking on the nucleons in the in- 
termediate states as well as the dressing of the D-meson, nucleons and pions. 



The effective DN interaction in the nuclear medium or G-matrix is obtained from the 
bare DN interaction, in a similiar way as done in Ref. [35] for the KN interaction. However, 
in this case, only the in-medium s-wave DN interaction is considered taking, as a bare 
interaction, a separable potential model as done for the K~p interaction in Ref. [36]. For 
the separable potential we use the following ansatz in momentum space 



where g is the coupling constant and A the cutoff. These two parameters will be determined 
by fixing the position and the width of the Ac(2593) resonance, the analogous to the A(1405) 
resonance in the charm sector. For the interaction matrix Qj, we use the standard result 
derived from SU(3) fiavor symmetry ( [36,37]). This bare interaction allows for the transition 
from DN to other channels, namely, ttAc, ttSc, jyAc and jyEc, all having charm c — 1. 
Therefore, we are confronted with a coupled-channel problem. In principle, one should use 
SU(4) symmetry as we incorporate charmed mesons. However, channels with a strange and 
charm quarks such as K'E^i DgA or DgT, have not been considered due to the fact that their 
respective masses are well above the DN threshold (see Fig. 1). Therefore, we only consider 
channels with up, down and charm-quark content keeping the SU(3) symmetry. 

The resultant meson-baryon G-matrices can be grouped in a matrix notation, where each 
box corresponds to one channel. The DN channel can have isospin J = or / = 1. In the 
first case, it couples to the ttSc and rjAc channels and the corresponding matrix has the 
following structure 



A. In-medium s-wave DN interaction 



Vij{k,k')^g^QjVi{k)vj{k') 

= |^G„e(A-fc)e(A-fc'), 



(1) 




while for / = 1 it 



can couple to the ttAc, ttEc and rjHc channels 




■DN ^■kKo-*DN <J-7rSc 
■ttAc G^ttAc— »7rAc Gt^y.c- 

•?7Sc GttA^^jjJ:^ G^rEc 
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Keeping this structure in mind, the G-matrix is formally given by 

(Mi^i I G{n) I M2B2)= {MiBi I V I M2B2} 

+ {MiB, I V I MsBs)——^^^^^§^—{M^B^ \ G{Q) \ M^B^) . (2) 

M3B3 " ~ - J^Bs + 

In Eq. (2), Mj and Bi represent the possible mesons {D, n, rj) and baryons (A^, Ac, Ec) 
respectively, and their corresponding quantum numbers such as spin, isospin, charm, and 
linear momentum. The function Qms.Bs stands for the Pauli operator which allows only 
intermediate nucleon states compatible with the Pauli principle. The energy variable is 
the so-called starting energy while ^/s is the invariant center-of-mass energy, i.e., ^/s = 
v/TF^^, where (Q, P ) is the total meson-baryon four momentum in a frame in which 
nuclear matter is at rest. 

The former equation for the G-matrix has to be considered together with a prescription 
for the single-particle energies of all the mesons and baryons participating in the reaction 
and in the intermediate states. These energies can be written as 

EMi{B,){k) = ^A;2 + m2^.(5.) + UM^{Bi){k, E^.^^,^) , (3) 

where UMi{Bi) is the single-particle potential of each meson (baryon) calculated at the real 
quasiparticle energy -E^.^^.^. For baryons, this quasiparticle-energy is given by 

EZ (k) = + m% + Re Us, {K E^) , (4) 
while, for mesons, it is obtained by solving the following equation 

W,(fc)f = fc' + mL. + RenM,(fc,^IfJ . (5) 

In the present paper we have considered the single-particle potential for the D-meson, 
nucleons and pions together with the decay width of the meson. The reason for the 
inclusion of this width will be clarified in the results section. 

For nucleons, as done in Ref . [38] , we have used a relativistic a — u model, where the 
scalar and vector coupling constants, Qs and respectively, are density dependent [39] . 

One of the most important modifications comes from the introduction of the pion self- 
energy, n^(/i;,cj), in the intermediate TrA^ and vrSc states present in the construction of the 
effective DN interaction. The pion is dressed with the momentum and energy-dependent 
self-energy [40] which was also studied in Ref. [38] . This self-energy incorporates a p-wave 
piece built up from the coupling to Ip — Ih, 1 A — Ih and 2p — 2h excitations, together with 
short-range correlations. The model also contains a small and constant s-wave part. 

The D-meson single-particle potential in the Brueckner-Hartree-Fock approach is 
schematically given by 

Unik, Ef) ^^{DN\ GoN^DNin = E'^ + Ef) \ DN), (6) 

N<F 

where the summation over nucleon states is limited by the nucleon Fermi momentum. As 
it can be easily seen from Eq. (6), since the effective DN interaction (G-matrix) depends 
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on the D-mcson single-particle energy, which in turn depends on the D-meson potential, 
we are confronted with a self-consistent problem. We also note that the G-matrix in the 
above equation becomes complex due to the possibility of DN state decaying into the ttAc 
and ttEc channels. As a consequence, the potential C/i? is also a complex quantity. Further 
details are given in the next section. 

The G-matrix equation for the s-wave in-medium DN interaction in the partial wave 
basis using the quantum numbers of the relative and center-of-mass motion reads 



{{MrBr)-k"\G\P,Q)\{M2B2);k)= {{MrB^)-k"\V'\{M2B2)]k) 
E l^dk'{{M,B,y,k"\V'\{MsBsy,k') 



+ 

M3B3 ■ 



x{{M3Bsyk'\G'{p,n)\{M2B2yk) , (7) 

where the variables k, k', k" arc the relative momenta and P is the linear center-of-mass 
momentum. The functions A;^ and k\^ are, respectively, the square of the momentum of the 
baryon and that of the meson in the intermediate states, averaged over the angle between 
the total momentum P and the relative momentum k' (see appendix A in Ref. [35]) 



kl{k',P)^k''+( 7 )\\ 



kl{k',P)^k''+( '7 Yp^ (8) 
VmM + msJ 

The angle average of the Pauli operator, Q M-isSk' i P) i differs from unity only for the DN 
channel (for details see appendix A in Ref. [35]). 



B. Z)-meson single-particle energy in the Brueckner-Hartree-Fock approximation 

For the s-wave component of the DN interaction, the Brueckner-Hartree-Fock approxi- 
mation to the single-particle potential of a D-meson embedded in a Fermi sea of nucleons 
[Eq. (6)] becomes 

UoikD, E'Eiko)) = ^ E(2^ + + 0' I k'dk f{k, ko) 

x{{DNyk\G'(P^,E'j^{kn)+E]^{M;)y{DNyk) , (9) 

where and k% are the square of the center-of-mass momentum and nucleon momentum, 
respectively, averaged over the angle between the external D-meson momentum in the lab 
system, kiy, and the DN relative momentum, k, used as integration variable in Eq. (9) (see 
appendix B in Ref. [35]). These angle averages eliminate the angular dependence of the 
G-matrix and allow to perform the angular integration in Eq. (9) analytically, giving rise 
to the weight function, f{k,kD), 
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f{k,k 



D 



1 for A; < ^rf^, 



ior \^kD - {I + Ok\ > kp, (.f.^ 

4^(1 + O^i^^ 

where ^ = — — and /ci? is the Fermi momentum. The magnitude of the relative momentum 
mo 

k is constrained by 

,_ = ^£_±|£. (U) 

After self-consistency for the on-shell value UoikD, E]^) is achieved, one can obtain the 
complete energy dependence of the self-energy Ilj:){kD,i-o), 

Uoiko, ou) = 2^kl + ml Uoiko, uj) (12) 

by replacing in Eq. (9) by u. This self-energy can then be used to determine the 
D-meson single-particle propagator in the medium, 

DD{kD,uj)^ ] , (13) 

(j'^ - kl-mjj- 2^ ml + kyjoiko^oj) 

and the corresponding spectral density 

Soiko^^^) = --l-mDoikD.^^) ■ (14) 

TT 

In fact, only the value of the potential L^d at the quasiparticle energy u = E]^ is de- 
termined self-consistently. This amounts to taking, in the subsequent iterations leading to 
self-consistency, the so-called "quasiparticle" approximation to the D-meson propagator 

Dmn^u^) = i— — , (15) 



kl-ml- 2Jml + klUo^kn, E'E) 



which gives rise to a simplified spectral strength 

S'E{kD,Uj) = 



2^ml + kllmUDikD,E'g) 

(16) 



^\u^-kl-ml- 2^ ml + klRe UoikD, E'E) P + | 2^ ml + kllm UoikD, E^) P 

The location and width of the peak in this distribution are determined, respectively, by the 
real and imaginary parts of UoikD, E^)- 

This self-consistent scheme was also used in previous works for the K meson [35,38] 
and, although it represents a simplification with respect to the more sophiscated scheme 
followed in Refs. [41,42] for the K case where the full energy dependence of the K self- 
energy is self-consistent determined, the approximation is sufficiently good as already shown 
in Refs. [35,38]. 
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III. RESULTS 



We start this section by showing the mass distribution of the ttSc state in Fig. 2. The 
mass distribution is given by 

^^C|T,«.,J^,,„ (17) 

where C is related to the particular reaction generating the vrSc state prior to final state 
interactions, pcM is the ttSc relative momentum and T^s°-^7rSc ^he isospin zero component 
of the on-shell s-wave T-matrix for the ttEc channel. The study of the mass distribution of 
the ttEc state in the 7 = channel reflects the Ac(2593) resonance as seen in Fig. 2. In this 
figure the mass distribiition of the ttSc state for J = is displayed as a function of the CM. 
energy for different sets of coupling constants g and cutoffs A. Wc observe that our coupled- 
channel calculation generates dynamically the Ac(2593) resonance. The position (2593.9 ±2 
MeV) and width (F = 3.6^1;° MeV) are obtained for a given set of couphng constants and 
cutoffs in the range between 0.8 and 1.4 GeV. This resonance was also obtained in the 
framework of the x-BS(3) approach based on the chiral SU(3) Lagrangian and formulated in 
terms of the Bethe-Salpeter equation [43]. It is interesting to observe that the set of coupling 
constants and cuttoffs that generate the resonance are very similar to the ones used in the 
s = — 1 sector to reproduce dynamically the A(1405) resonance [36]. 

Once the position and width of the Ac(2593) resonance are reproduced dynamically, 
we study the effect of the different medium modifications on the resonance. The real and 
imaginary parts of the resulting in-medium s-wave DN amplitudes in the 1 = channel for 
A = 1.0 GeV and gi^=13.4 and for a total momentum \ kn + kN \= are given in Fig. 3 as a 
function of the invariant center-of-mass energy for different approaches: T-matrix calculation 
(dotted lines), including Pauli blocking (dot-dashed lines) and self-consistent calculation for 
the D-meson (solid lines) at nuclear matter saturation density po = 0.17 fm~^. Due to the 
fact that the resonance lies just few MeV above the vrSc threshold, we have included the 
decay width of the Eg baryon in our calculations by considering the corresponding imaginary 
part for the Eg self-energy. In this way, the resonance is reproduced in a more realistic way, 
even below the ttEc state. We clearly see, as noticed already for the A(1405) resonance, the 
repulsive effect of the Pauli blocking on the resonance, being generated at higher energies. 
However, as already pointed out for the KN interaction, the shift in energy of the Ac(2593) 
resonance is intimately connected with the energy dependence of the DN interaction and 
changes are expected from a self-consistent incorporation of the Z^-meson properties in the 
DN interaction. Therefore, the effect of dressing the D-meson with the complete complex 
self-energy on the in-medium DN amplitude makes the resonance peak stay pretty close or 
even at lower energies than its free space location. This is due to the attraction felt by the 
D-mcson that compensates the repulsive effect induced by Pauli blocking on the miclcon, as 
already seen for the A(1405) resonance in Refs. [35,36,41,42,44]. On the other hand, the real 
and imaginary parts of the DN amplitude become much smoother because of the £)-meson 
strength being spread out over energies. 

In Fig. 4 we display the previous plot of the real and imaginary parts of the s-wave 
/ = DN amplitude together with the 1 = 1 channel for a larger energy scale and a smaller 
range in the y-axis. Apart from the previous three approaches, T-matrix calculation (dotted 
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lines), in- medium calculation only including Pauli blocking effects (dot-dashed lines) and 
self-consistent calculation for the D- meson potential (solid lines), we also include the in- 
medium properties of the nucleons together with the pion self-energy in the self-consistent 
process for the D- meson potential (long-dashed hnes) . While the nucleon dressing only shifts 
the Ac(2593) resonance in energy, as seen for the A(1405) resonance, the pion dressing could 
introduce important changes to the previous self-consistent procedure. In this approach 
and for 7 = 0, the imaginary part of the DN interaction becomes smoother in the region 
of energies where the Ac(2593) resonance was generated in the previous three approaches. 
For energies below the ttEc threshold, we observe a bump around 2.5 GeV that, in our self- 
consistent many-body approach, can decay to states such as 7r(p/i)Ec or D{h.Jn:)N , where 
in parentheses we have denoted an example for the component of the vr and D-meson that 
show up at energies below ttSc- On the other hand, it is also interesting to observe a second 
structure in the 7 = amplitude, already seen for the three initial approaches, which lies 
below the DN threshold of 2.806 GeV. This feature in the imaginary part around 2.75 GeV 
seems to be enhanced in the full self-consistent procedure. This structure, like the previous 
one, is also a state with the Ac-like quantum numbers and, in this case, it can decay to the 
ttSc state. Whether the first resonant structure is the in-medium Ac(2593) resonance and 
the second bump is a new resonance is something that deserves further studies. Note that 
the energy region of interest for the calculation of the D-meson potential for / = lies on 
the right-hand side of the free Ac(2593) resonance where this second structure shows up for 
all approaches. Actually, this is a different feature of the DN interaction compared to the 
case of the KN interaction. The KN interaction was basically determined by the behaviour 
of the s-wave / =0 A(1405) resonance in the medium. For the 1 = 1 component of the 
DN interaction, the self-consistent calculation tends to dilute any structure present in the 
T-matrix. It is seen that for the full self-consistent calculation, the 7 = 1 component is 
considerably reduced. A more detailed analysis concerning the DN interaction and, hence, 
the Z)-meson potential is presented in the following plots. 

In Fig. 5 we plot the real and imaginary parts of the D-meson potential at /cd = as a 
function of density for A = 1 GeV and =13.4 and for the three previous in- medium cal- 
culations: in- medium calculation with only Pauli blocking (dot-dashed lines), self-consistent 
calculation for the D-meson (solid lines) and self-consistent calculation for the D-meson 
including the dressing of nucleons and the self-energy of pions (long-dashed lines). The D- 
meson potential becomes more attractive as we increase the density in all three approaches, 
showing a less smooth behaviour in the case when nucleons and pions are dressed both in 
the real and imaginary parts. However, the potential turns out to be more attractive for a 
self-consistent calculation with respect to the case when only Pauli blocking is considered. 
If we compare these results to the case of the K potential, the picture depicted here is 
different. For the K meson case, the shift of the A(1405) resonance in energy due to the 
Pauli blocking changes the scattering amplitude at the threshold from being repulsive 
to attractive giving rise to an attractive K potential. However, when the properties of K 
were incorporated self-consistently, the attraction was drastically reduced [35,41,42]. For 
the D-meson, the DN threshold is 2.806 GeV and, therefore, the DN amplitude is studied 
for energies on the right-hand side of the free Ac(2593) resonance, away from this resonant 
structure. Furthermore, in the K meson case, the s-wave 7 = amplitude turns out to 
be the main contribution to the K potential due to the presence of the A (1405) resonance. 
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For the D-meson, as the energy region of interest is not sitting on the resonance, the I — 1 
could also become important in the calculation of the D-meson potential. 

In order to study the isospin dependence of in-medium DN interaction in more detail 
as well as the dependence on the chosen set of coupling constants and cutoffs, in Fig. 6 and 
7 we represent the real and imaginary parts of the DN amplitude at p — po in the 7 = 
and I — 1 channels as functions of the center- of- mass energy around the DN threshold at 
total momentum | ko + k^ |= for the self-consistent calculation when only the D-meson is 
dressed (Fig. 6) and for the full self- consistent calculation (Fig. 7). Furthermore, in Fig. 8 
we show the real and imaginary parts of the D- meson potential at A;£) = as a function of 
the density for A = 1 GeV and g"^ =13.4 including the isospin decomposition. 

For the first self-consistent approach where only D-mcsons are dressed ai p — po (Fig. 6), 
the real parts of the / = and 1 = 1 components turn out to be of the same order of 
magnitude. Therefore, it is not clear which isospin will determine the behaviour of the D- 
meson potential. Actually, the potential depends on the exact value of both contributions 
taking into account the isospin factor 27-1-1 for each isospin component (see Eq. (9)). For 
A = 1 GeV, the main contribution for the potential at p = po comes from the J = 1 
component, as seen in Fig. 8. On the other hand, we note that for A = 0.8 GeV, the 1 = 
component moves from attraction to repulsion as the energy increases (Fig. 6). Nevertheless, 
the region of integration to obtain the £)-meson potential for all the sets of parameters of 
couphng constants and cutoffs lies around 2.80-2.84 GeV at p = po and, therefore, we obtain 
a repulsive effect for / = and an attractive effect for / = 1 in the potential at p = po for 
all sets (see Fig. 8 for A = 1 GeV). The imaginary part displays a smoother behaviour than 
in the case of the full self-consistent calculation. When the full calculation is performed 
(Fig. 7), the peak observed around energies of 2.75 GeV gets diluted with increasing cutoff. 
However, the real part shows more structure than in the previous self-consistent approach 
and the 7 = component becomes more important with respect to the previous case. This 
effect together with a less important contribution of the 7 = 1 component makes the 7 = 
govern the 7}-meson potential for densities larger than p = 1.5po ior all sets. For p = po 
and A = 1 GeV, the attractive 7 = 1 component still dominates. What it is not so clear is 
whether the 7 = will be attractive or repulsive for all the sets as the region of integration 
lies around 2.71-2.78, where the 7 = real part changes in sign. The different region of the 
integration in this last case with respect to the case where only the 7^-meson is dressed is 
mainly due to the inclusion of an attractive potential for the nucleons as the value of the 
7^-meson potential does not change by a great amount for both self-consistent processes (see 
Fig. 8). 

In Fig. 8 we observe that, in the case where only the 7^-meson is dressed and for A = 1 
GeV, the real part of the 7^-meson potential is governed by the 7 = 1 component as density 
grows. This is also seen for the other set of parameters studied. On the other hand, for the 
full calculation, the 7 = component controls the behaviour of the 7^-meson potential as we 
increase the density not only for A = 1 GeV but also for the other sets. This is due to the 
resonant structure seen in Fig. 7 that also causes a less smooth behaviour of the 7)-meson 
potential. 

In order to study the dependence on the cutoff and coupling constant, we show in Fig. 9 
the 75-meson potential at /cd = as a function of the density for different sets of coupling 
constants and cutoffs and for the two approaches mentioned before. When only the 7?-meson 
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is dressed, higher densities are required in order to obtain an attractive potential with respect 
to the full self-consistent calculation. In both cases, the D-meson potential tends to get more 
attractive as density increases. The reason is that, in order to obtain the D-meson potential 
as density increases, we are integrating over a larger energy region where the interaction is 
attractive. However, as mentioned in the previous plot, the isospin component responsible 
for that behaviour is different depending on the approach. On the other hand, the value 
of the potential turns out to be slightly sensitive to the chosen set of parameters. There 
is, however, in both cases a density for which all sets give the same value for the D-meson 
potential. When only the D-meson is dressed, we obtain a range of values for the D-meson 
potential at p = po between 8.6 MeV for A = 0.8 GeV and -11.2 MeV for A = 1.4 GeV. 
For the full self-consistent calculation, the range of values covered lies in between 2.6 MeV 
for A = 0.8 GeV and -12.3 MeV for A = 1.4 GeV. This result would indicate that our self- 
consistent coupled- channel calculation gives a different prediction for the D-meson potential 
as compared to the calculations based on the QGD sum-rule approach [16,19], the quark- 
meson couphng model [21], lattice calculations [20,31,32] and a recent work based on a chiral 
model that obtains the medium modification of the mass of D-meson due to its interaction 
with the light hadron sector [33]. The calculations based on the QGD sum rule approach 
as well as quark-meson coupling model predict the mass drop of the D-meson to be of the 
order of 50-60 MeV at nuclear matter density. A similar drop at finite temperatures [20] is 
suggested from the lattice calculations for heavy quark potentials. An interaction Lagrangian 
based on chiral perturbation theory also gives a similar shift at nuclear matter density, when 
the Tomozawa- Weinberg interaction is supplemented by the attractive DN Sigma term and 
the charm content of the nucleon is ignored [33]. On the other hand, the effective chiral 
Lagrangian approach [33] is seen to yield a larger drop of the D-meson masses (~ 200 MeV 
at nuclear matter density). In the present calculation, the coupled-channel effects seem to 
result in an overall drastic reduction of the in-medium effects independent of the parameters 
{g,A) and the in-medium properties of the intermediate states compared to the previous 
works. 

Finally, once self- consistency is reached, we calculate the full energy dependence of the 
£)-meson self-energy which defines the in-medium £)-meson single-particle propagator and 
its spectral density (see Eqs. (13-14)). The spectral density at zero momentum is shown in 
Fig. 10 for A = 1 GeV and for several densities in the two approaches considered before. In 
the case when only the D-meson is dressed in the self-consistent process, the peak of the 
D-meson pole moves towards lower energies as density increases since the D-meson potential 
becomes more attractive for A = 1 GeV (see Fig. 9). Moreover, it is observed that the D- 
meson spectral density falls off more slowly on the left-hand side of the quasiparticle peak 
as density increases. This is due to the 1 = 1 component of the off-shell D-meson self- 
energy, which is related to the D-meson potential through Eq. (12). As density increases, 
we have already discussed that the 1=1 component governs the behaviour of the potential 
and, hence, the self-energy. Therefore, any structure lying close the quasiparticle peak is 
enhanced and clearly visible in the spectral density. On the other hand, we also notice 
some structure of the spectral density to the left of the quasiparticle peak at energies of the 
D-meson of around 1.63-1.65 GeV, the origin of which could be traced back to the presence 
of the Ac(2593) resonance. 

For the full self-consistent calculation, the quasiparticle peak moves towards lower ener- 
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gies as density increases according to the behaviour of the D-meson potential seen in Fig. 9, 
although there is region of densities around 0.5po where the D-meson potential turns out 
to be repulsive. For densities between po and 1.5 Pq, the quasiparticle peak mixes with a 
structure, already noticed on the left-hand side of the quasiparticle peak for 0.5 po, mak- 
ing it more difficult to separate the contribution of the quasiparticle peak and the effect of 
this structure. This bump in the spectral density is due to the structure observed in the 
imaginary part of the / = component of the in-medium DN interaction which lies close 
to the DN threshold. The nature of this structure and in-medium changes of the D-meson 
production could probably be studied experimentally in the near future with the PANDA 
experiment at the GSI International Facility in Darmstadt [45] . 

IV. CONCLUSIONS 

We have performed a microscopic self-consistent coupled-channel calculation of the 
single-particle potential of the D-meson and, hence the D-meson spectral density, embedded 
in symmetric nuclear matter assuming a separable potential for the s-wave DN interaction. 

The Ac(2593) resonance, which is the counterpart of the A(1405) in the charm sector, 
has been generated dynamically for a given set of coupling constants g and cutoffs A. 

We have also studied the medium effects on the Ac(2593) resonance and, hence on the 
D-meson potential, due to the Pauli blocking of the intermediate nucleonic states as well 
as due to the dressing of nucleons and pions. The D-meson optical potential has been ob- 
tained for different approaches: in-medium calculation including only Pauli blocking on the 
intermediate nucleonic states, self-consistent calculation for the D-meson and self-consistent 
calculation for the D-meson including the dressing of nucleons and the pion self-energy. 

We observe that the two self-consistent schemes show a stronger density dependence 
compared to the case when only Pauli blocking effects are included. When only the D- 
meson is dressed in the self-consistent procedure, the D-meson potential at p = po stays 
between 8.6 MeV for A = 0.8 GeV and -11.2 MeV for A = 1.4 GeV. For the full self- 
consistent calculation, the range of values covered lies in between 2.6 MeV for A = 0.8 GeV 
and -12.3 MeV for A = 1.4 GeV. We conclude that the coupled-channel effects seem to result 
in an overall reduction of the in-medium effects independent of the set of parameters {g,A) 
and the in-medium properties of the intermediate states compared to previous work. 

The isospin dependence of the D-meson in both approaches have been also analyzed. 
While the 1 = 1 amplitude governs the behaviour of the D-meson potential when only D- 
mesons are dressed, the D-meson is basically determined by the / = component in the full 
self-consistent calculation. 

The D-meson spectral density has been finally obtained for both self-consistent schemes. 
Although the quasiparticle peak stays closer to its free position in both cases for nuclear 
matter saturation density, the features of the low-energy region on the left-hand side of 
the quasiparticle peak are different according to the different in-medium behaviour of the 
Ac(2593) resonance in both approaches. 

The in-medium effects devised in this work can be studied in heavy-ion experiments at 
the future International Facihty at GSI. The PANDA experiment at GSI [45] will measure 
hadrons with charm by antiproton beams on nuclei with its microvertex detector. In-medium 
changes of open charm hadrons can be addressed by the study of the excitation function 
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and the correlation function of D- and D-mesons. We stress that, ahhough the in-medium 
potential for D-mesons has turned out to be quite small in our investigation, the production 
of D-meson in the nuclear medium will be still enhanced due to the additional strength of 
the D-meson spectral function below the quasiparticlepeak. This effect is similar to the one 
extracted for the enhanced K production in heavy- ion collisions [46]. 

The present coupled- channel approach to the D-meson properties in the nuclear medium 
is, as the first of its kind, exploratory and can be improved by incorporating relativistic 
potentials as well as chiral constraints on the bare hadronic interactions. It will be also in- 
teresting to explore the in-medium effects for D-mesons in dynamical approaches for studying 
e.g. the excitation function, which we will leave for future work. 
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